Logarithmic bred vectors in spatiotemporal chaos: structure and growth 
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Bred vectors are a type of finite perturbation used in prediction studies of atmosplieric models 
that exhibit spatially extended chaos. We study the structure, spatial correlations, and the growth- 
rates of logarithmic bred vectors (which are constructed by using a given norm). We find that, 
after a suitable transformation, logarithmic bred vectors are roughly piecewise copies of the leading 
Lyapunov vector. This fact allows us to deduce a scaling law for the bred vector growth rate as 
a function of their amplitude. In addition, we relate growth rates with the spectrum of Lyapunov 
exponents corresponding to the most expanding directions. We illustrate our results with simulations 
of the Lorenz '96 model. 
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I. INTRODUCTION 

In dynamical systems theory, chaos is usually defined 
on the basis of the exponential departure of infinitesi- 
mally separated initial conditions. The exponential law 
stems from the linearity of the equations that govern in- 
finitesimal perturbations. In applications, however, er- 
rors are typically finite. This occurs, for instance, in the 
breeding method used to generate finite perturbations 
for ensemble forecasting at the National Centers for en- 
vironmental Prediction (USA) [l|, 0]. Bred vectors are 
closely related to Lyapunov vectors (LVs) albeit different 
since they are finite by construction and result from the 
evolution of perturbations, which are imposed to have a 
certain size via periodic normalizations. One advantage 
of bred vectors in applications is that the model under 
study does not have to be linearized, in contrast with 
LVs. (Keep in mind that the linearization of a meteoro- 
logical model is a delicate question due to the presence of 
nondifferentiable or discrete programming structures, see 
e.g. [3|-) The drawback is that bred vectors are governed 
by fully nonlinear models, what constitutes a challenge 
for their theoretical description. 

In this article we use a special class of bred vector, the 
so-called logarithmic bred vector (Log-BV) [J, |5[ , which 
has the particularity of being normalized using the ge- 
ometric (or zero) norm. This choice has proven to be 
the most convenient for spatially extended systems (see 
Ref. [3, Q) and offers some advantages in theoretical 
terms. We will focus on two aspects of bred vectors: 
First, we uncover common structural properties of Log- 
BVs and the leading LV. In the second part of the paper 
we study the growth rate of Log-BVs, which we shall 
refer to as bred exponents (BEs) — instead of the more 



imprecise term finite-size Lyapunov exponents [6| — for 
the sake of brevity and to emphasize their origin. We 
find an interesting relation between BEs and Lyapunov 
exponents (LEs) of spatio-temporally chaotic systems. 

This paper is organized as follows: In Sec. |TT]we intro- 
duce the Lorenz '96 model, which we will use through- 
out this contribution in order to illustrate our consid- 
erations. For the parameter values selected the Lorenz 
'96 is hyperchaotic and exhibits spatio-temporal chaos. 
In the following section we explain the computation of 
Log-BVs. A surface picture is introduced for the lead- 
ing LV in Sec. IIV Al and extended to analyze Log-BVs in 
Sec. lIVBl In Sec.|V]we study the growth rates [i.e., BEs) 
of Log-BVs. In Sec. IV Al the convergence to the largest 
LE is found to obey a scaling relation. In Sec. IV Bl we 
explain how the BEs are related to the spectrum of LEs. 
The conclusions are summarized in Sec. I VII 



II. THE LORENZ '96 MODEL 

The Lorenz '96 model [3, Q is a toy model originally 
proposed in the context of atmospheric dynamics and 
used extensively to test novel techniques and applica- 
tions. It is a time continuous model consisting of a set of 
nonlinear ODEs coupled in a ring geometry: 



du{x,t) 
Jt 



—u{x, t) — u{x — 1, t)u{x — 2, t) 

+u{x - 1, t)u{x + 1, t) + F, 
with X — \, ..., L. (1) 
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u{x,t) can be seen as a scalar meteorological variable, 
e.g. temperature, at L equally spaced sites x on a lati- 
tude circle (and hence periodic boundary conditions are 
assumed). Moreover, there is an external forcing con- 
stant F that mimics the solar driving of the atmosphere. 
For i > 3 the solutions of ([1]) are chaotic if F is large 
enough. In particular the steady solution u{x,t) = F 
becomes unstable if F > 8/9 [8] . A more detailed study 
shows that stable nonchaotic solutions survive up to a 



value of F that, though depending on L, is approximately 
in the range 4 to 6 ;9'] . Beyond some threshold value of F, 
chaotic dynamics of the model becomes fully developed. 
More precisely, the dynamics is extensive with L (see 
e.g. results in [lOj for F — 10). Extensivity means that 
many relevant quantities (dimension, entropy, etc.) scale 
linearly with the system size, and the LEs {Xn}n=i,....L 
converge to a density in the "thermodynamic" limit L ^■ 
oo. This property is shared by a number of extended 
dynamical systems ranging from coupled map lattices to 
PDEs such as the Kuramoto-Sivahinsky or the complex 
Ginzburg-Landau equations [llj . 

In order to compute bred and Lyapunov vectors and 
their growth rates in the Lorenz '96 model, we integrate 
Eq. ([1]) and its linearization (tangent space) by using a 
fourth order Runge-Kutta-solver with time step At = 
0.01. Before measuring the quantities we are interested 
in, we allow the system to breed Log-BVs for a transient 
time ttrans > 8(L^-^). 



III. LOGARITHMIC BRED VECTORS AND 
BREEDING 

In 1993, Toth and Kalnay [l[ created a special oper- 
ational cycle designed to "breed" fast growing errors in 
meteorological models. It uses finite perturbations that 
are periodically normalized to become bred vectors after 
some breeding cycles. 

In the following we describe the procedure to compute 
bred vectors in mathematical terms. A control (unper- 
turbed) trajectory u(t) = [u{x,t)]^'Zi and a perturbed 
one u(t) = [u{x, t)]^Zi are evolved in parallel obeying 
Eq. ([T]). The difference between the perturbed and un- 
perturbed trajectories is calculated every time interval r, 
say at times Tm = mr, with m = 1, 2, ... to obtain 



6u{t„i) = u(t,„) - u(r„ 



(2) 



These differences are then scaled down at i = r^ to a 
given perturbation amplitude Eq by defining 



SUsoiTrn) = So 



IIMt™)! 



(3) 



where Sui;g{x,Tm) denotes the bred vector at time t,„, 
and is || • || a particular norm. The perturbed trajectory 
is then redefined by means of the bred vector: 



"('^m) = u(Tm) + 5Ueo(T,„), 



(4) 



with T+ = liniix^o Tm + i^ referring to the same time r^, 
but after the rescaling. Perturbed and unperturbed tra- 
jectories are again integrated forward in parallel until 
next rescaling scheduled at time Tm+i- 

The definition of the bred vector contains two ingredi- 
ents. One is the parameter eg controlling the amplitude. 
The second ingredient is the norm to be used since differ- 
ent norms will produce different bred vectors. In recent 



works [3, ^ it was found that the 0-norm (or geometric 
norm) 



\\Su{t)\\ ^ hm 

9— >0 



nl/« 



J2\Suix,t)\'' 



x = \ 



= ^\8u{x,t)\ 



l/L 



x=l 



(5) 

is a convenient choice for breeding. Due to the multiplica- 
tive character of the linear dynamics the 0-norm ought to 
produce bred vectors that, at different times, Tm are the 
most statistically equivalent among them. The bred vec- 
tors constructed in this way are called logarithmic bred 
vectors (log-BVs) 0,11 d. 



IV. STRUCTURE OF LOGARITHMIC BRED 
VECTORS 

A. The main Lyapunov vector 

Before starting our analysis of the Log-BVs it is use- 
ful to recall recent results concerning the structure of the 
main (or leading) LV. In a dynamical system an infinitesi- 
mal perturbation evolves generically towards the leading 
Lyapunov vector g(t), which indicates the direction of 
maximal growth for perturbations integrated since the 
remote past. The orientation of the LV in tangent space 
depends on the position in the chaotic attractor (that can 
be parametrized by time). For extended chaotic systems, 
it is observed that the LV projects very inhomogeneously 
on space. More precisely, the vector localizes at some 
quite narrow region of the system [l3[. But noticeably, 
this localization is dynamic and the localization center 
changes as time evolves, thus recovering the homogeneity 
of the system in a statistical sense. In the 90s Pikovsky 
and coworkers [ij, \l^ found very useful for the theoreti- 
cal analysis to associate a "surface" with the leading LV 
by means of a Hopf-Cole transformation, which in d = 1 
dimension reads: 

hLvix,t) ^ln\g{x,t)\, (6) 

with g(i) — [g{x,t)]%'Zi. For a large family of systems 
p^ . including the Lorenz '96 model [Iq], hLv{.x,t) ex- 
hibits correlations in space and time which are described 
by the canonical Kardar-Parisi-Zhang (KPZ) equation of 
stochastic surface growth [I^. This mapping leads to 
interesting scaling properties for hLv{x^t). The average 
width w = {(hhy — (ft-Lv))^)^^^ scales with the length 
of the system [l^, [l^ as w ex L"; with a = 1/2 as in 
the KPZ equation in one dimension. This means that, at 
sufficiently long scales, /i^y (x, t) appears as the path of a 
random walk in d = 1. This self-affine profile translates 
into a power-law dependence of the structure factor (a 
spatial power spectrum) at small wavenumbers: 



S{k) oc fc-(2"+'*) for k<k, 



(7) 



where S{k) = {h(k,t)h{—k,t))t, and h{k,t) = 
^^ exp(jfca;)/i(a:,i). At short length scales, k > k, non- 
universal short-range correlations are expected to appear 
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FIG. 1: (color online) Upper panel: Surfaces associated with 
the LV and with the Log-BV with the smallest value of eo 



used e™ 



(Lorenz '96 model with F = 8.0 and L 



1024). Lower panel: Snapshot of the difference fields »7e„ = 
h^g {x, t) — hLv{x, t) (curves are arbitrarily shifted to improve 
their visibility). 



due to the deterministic character of the system. How- 
ever, below k universal scaling properties emerge and this 
is reflected in a spatial correlation that generically decays 
as ~ k~^ in one dimension. 



B. The logarithmic bred vectors 

Analogously to the surface associated with the leading 
LV in Eq. ([6]), one can define a surface 



heo{x,t) ^\n\5u^g{x,t)\ 



(8) 



by Hopf-Cole transforming the Log-BV of amplitude Eq. 
The main LV evolves following linear equations and, 
therefore, only its direction in tangent space matters. 
However for a bred vector the norm plays a very impor- 
tant role through the periodic rescaling of the perturba- 
tions defined above. If the geometric norm is chosen, the 
spatial mean of the associated Log-BV surface is fixed to 



- 1 

h= —^h^g{x,t^ Tm) = Ineo- 



(9) 



x = l 



In the limit of small amplitude Eq -^ 0, the Log-BV aligns 
with the leading Lyapunov vector, and their profiles co- 
incide, h^g^Q{x,t) — hLvix,t), apart from an arbitrary 
constant due to the arbitrary norm (i.e. arbitrary h) of 
the leading LV. This fact suggests that the leading LV 
can be a good reference point to analyze the structure of 
Log-BVs. In Fig. [1] we show a snapshot of the difference 
fields 



Veo{x,t) = Koi.^,t) -hLv{x,t). 



(10) 



We can see that, unless Eq is too large, rj^^ is mainly 
composed of plateaus, which indicates that the structure 
of h^g{x,t) is roughly a piecewise copy of hLv{x,t). Re- 
markably, this structure has been previously observed for 
LV-surfaces associated with LEs A„ {n > 2) smaller than 
the largest one A = Ai [H, HO] ■ For the Log-BVs the typ- 
ical plateau size decreases as eo is increased, whereas for 
the sub-dominant LVs the plateau size is known [Ig, [20| 
to decrease as the index n is increased (A„ > A„+i). This 
suggests an interesting relation between both (Lyapunov 
and bred) vector types that is exploited below in Sec. lVBI 
Note also that Primo et al. [5| used the surface ([8]) 
to uncover several spatio-temporal scaling relations that 
already revealed in an indirect way the existence of cut- 
off lengths. These cut-off lengths are now apparent in 
the light of the plateaus in Fig. [TJ 



V. BRED EXPONENTS 

We now focus on the (exponential) growth rate of Log- 
BVs, which we will refer to as bred exponents (BEs). We 
define the BEs as: 



A. 



_ 1 / \\du{T„^ + 
r\ ||(5u£(,(t„ 



(11) 



with r denoting the time between rescalings. X^g can 
be seen as a type of finite-size Lyapunov exponent, see 
Appendix A of [6] . The value of A^,, is not very sensitive 
to r if T remains small. In our simulations we chose r 
such that the perturbation does not amplify more than 
exp(l/2) times in a breeding cycle, so we take 



r^<4. 



(12) 



with A being the largest LE. First of all, it is convenient 
to transform Eq. (fTTj) . using Ec^. (|9]): 



A. 



[{h{x,Tm +r)) -Ineo] 



(13) 



This expresses in mathematical terms that the BE is the 
average velocity of the Log-BV-surface. In the transfor- 
mation of PT|) into ((T5)) we are assuming the geometric 



norm. Again, this choice makes plenty of sense because 
the geometric norm actually yields the least-fluctuating 
LE in systems with spatiotemporal chaos |15l |. 



A. Convergence of the Bred Exponent to the first 
Lyapunov exponent 

In this subsection we explore the dependence of the 
BEs, Aeo(L), on the amplitude £o and on the system size 
L. Note that in the limit of vanishing amplitude one 
recovers the largest LE of the system: Xi;^^q{L) = X{L) 

The key step of our following analysis is the use of the 
associated surfaces and the universal scaling laws they 
obey. Let us denote by X{L = oo) the LE of the model 
with infinite size. As already reported in Ref. [lal, the 
LE of a system of size L, X{L) deviates from the infinite 
size limit as 



X{L = oo) - X{L) 



1 

L 



(14) 



This stems from the fact that X{L) is the velocity of 
the associated surface, that scales as solutions of the 
KPZ equation: The asymptotic velocity of a KPZ- 
surface presents a system-size correction 2l| of order 
r^ 7^-2(i-q)^ with a = 1/2 one dimension. 

For the dependence of the BE on eo we resort to very 
simple arguments. Following the reasonings in [J| let us 
assume a Log-BV is (locally) linear whenever \Su{x, t)\ < 
B, for a certain bound B. This bound defines borders 
of regions of size Ic where the Log-BV-surface is approx- 
imately a copy of the leading LV-surface. The KPZ sur- 
face is self affine, and recalling Eq. ©, we expect the 
scaling law 



InS-lneo-Z" 



(15) 



to be fulfilled. Taking into account that the velocity is 
shifted from the asymptotic value by the inverse of the 
size, i.e. l^^: 



X{L = cxd) — Ae 



(InS-lneo) 



\/a 



(16) 



This relation is not the asymptotic one because when eo 
becomes extremely small the "nonlinear barrier" at In B 
is seldom achieved and the dominant correction is given 
by Eq. dHD. 

Scaling relations ((T4)) and (fT6)) can be cast into 



(17) 



X{L = ^)-X,,{L)^f{p)/L 



where p — (In S — lneo)^/i and /(p) is a scaling function: 
For p >• 1, f{p) = const., as expected from (fT4|) : whereas 
for p ^ 1, f{p) should scale as '^ p~^, according to ([TC|. 
Using the data from our simulations we can test the va- 
lidity of this scaling law. Figure [2] shows a very good 
collapse with just one fitting parameter, Ini? = 2.2, for 
different system sizes and values of Eq- Our scaling ar- 
gument is in good agreement with numerical data shown 
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FIG. 2: (color on line) Data collapse of X^q (L), defined in (|11|) 
and p3|l , for different values of £o and system size L {F = 8). 
Points fall approximately on a line, the scaling function g in 
(|17p . Inset: The value of the LE X{L = oo) was determined 
by extrapolating the largest LE from finite systems using (|14|) 



in Fig. [21 where one can see how the p ^ 1 asymptote 
converges to 1/p (shaded region in Fig. [5]). However, the 
approximation appears to be too crude to describe the 
behavior very close to the crossover point p w 1, where 
the functional form of /(p) deviates from the asymptote. 
Nonetheless the correct scaling variable p for the scal- 
ing relation ([T7]) and the two asymptotes are correctly 
captured. The crossover of /(p) at p « 1 marks the de- 
parture from a regime dominated by the finite-size sys- 
tem effects to the dominance of the amplitude Eq of the 
Log-BV. 



B. Bred exponents vs. Lyapunov spectrum 

In the preceding subsection we have introduced the 
concept of bred vector exponent X^g as the growth rate 
for a Log-BV of amplitude Sq. We have also seen how 
this BE approaches the leading LE in a finite size system 
as the amplitude eg is varied. Now we devote the present 
subsection to study the connection of the BEs with the 
spectrum of LEs. The question we want to address is 
to what extent growth rates of Log-BVs (finite errors) 
approach growth rates of LVs (infinitesimal errors). 



1. Main hypothesis 

From the scaling arguments discussed above, it can be 
expected that the BE X^o ought to be close to the n-th 
Lyapunov exponent A„ whenever both of them are piece- 
wise copies of the leading Lyapunov vector with similar 
plateau sizes. This hypothesis reads more formally: 



Aeo«A„ ff ?f^(eo)«^c (n) 



(18) 
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FIG. 3: (color online) The lines with symbols are structure 
factors Sik) obtained for Log-BVs with different perturba- 
tion amplitudes for F = 8.0 and L — 1024. The straight lines 
without symbols indicate the method of estimating the values 
fcc(eo) which are determined from the intersection of a hori- 
zontal line at S{kmin = 2'k/L) and a power-law S{k) oc fc~^"^ 
obtained from fitting the structure factor for rnin(lneo). 



where /^^ (Ic^) indicates the typical length scale over 
which the Log-BV (the n-th LV) surface is a piecewise 
copy of the leading LV. In the following we will drop the 
superindices BV and LV, as it is clear from the argument 
(eo or n) what vector type we are referring to. 



2. Estimation of crossover wavenumbers kc{eo) 



In Sec. IIVBI we showed that Log-BVs are piecewise 
copies of the main LV. This observation is now exploited 
to obtain a connection between BEs and LEs. In order 
to do so, we find it very convenient to look at the form 
of spatial correlations in Fourier space. The structure 
factors of Log-BV-surfaces with different values of eg are 
presented in Fig. [31 As we have seen in Fig. [1] the Log- 
BV-surface follows the leading LV-surface at short scales 
and, as a consequence, their respective structure factors 
should overlap above a certain fcc(£o)- This implies that 
both (Log-BV and LV) surface types share an interval 
kc{£o) < fc < fc — 0.3 with power-law structure factor 
S{k) r^ k^^. Log-BVs exhibit fiat structure factors for k 
smaller than a certain kdso), which indicates that dis- 
tant regions (corresponding to small k) are basically un- 
correlated. The crossover wavenumber fcc(£o) '~ ^c(eo)~^ 
is monotonically increasing with Eq and, from Eq. (1151) . 
we have fcc(eo) ^ {\nB - Ineo)"^/" with a == 1/2. No- 
tably, surfaces associated with LVs for n > 1 have a 
structure factor that also exhibits a knee at a certain 
crossover wavenumber kc{n) [16|, |20|. In sum, structure 
factors of LV-surfaces and Log-BV-surfaces are very simi- 
lar (cf. Fig. 4 in ^^), and they only differ in the algebraic 
dependence below the particular crossover wavenumber 
kc- 1/fc-type for the LVs, and flat (fc°) for the Log-BVs 



FIG. 4: (color online) Collapse of the structure factors using 
the the values fcc(£o) determined from the structure factors 
S{k) in Fig. El 
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FIG. 5: The lines with symbols indicate structure factors S{k) 
obtained for bred vectors with different perturbation ampli- 
tudes with F = 6,8, 14, 20 and L = 1024. The values fcc(eo) 
are determined like in Fig. [S] 



|12| , reflecting their different spatial correlations at large 
scales. 

The crossover wavenumbers fcc(£o) of the Log-BV- 
surfaces hgg{x,t) can be systematically extracted from 
their structure factors S{k), see Fig. |3] for L = 1024 
and F = 8. Figure H] shows the collapse of different 
structure factors when normalizing the wavenumbers by 
kc{£o), giving an idea of the goodness of the estimation 
of fcc(eo)- Figures [S] and [7] show that the same procedure 
can be used for other values of F and smaller systems. 



3. Numerical results 

In extended systems, it is customary to represent LEs 
with the index n normalized by the system size (or by the 
number of degrees of freedom) . This is done so to high- 
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FIG. 6: (color online) Lyapunov spectra (solid black lines) 
and bred exponents (symbols). Size L = 1024, and F = 
6,7,. ..,20. 



light the extensivity (assuming it exists) of the model 
because LEs for different systems sizes with identical pa- 
rameter values will approximately fall on the same line. 
Figure [H] shows the LEs connected by black lines. The 
Lyapunov spectra were obtained using the standard al- 
gorithm by Benettin et al. [22^i3] ■ 

We found in Ref. [20] for an extensive coupled map lat- 
tice that kc{n) ~ {n/ L) , with around one. This means 
that representing the LEs versus [n/L) is, for moderate 
values of n, equivalent to use a quantity proportional to 
kc{n) in the x-axis. 

If now we want to test our hypothesis that BEs and LEs 
are similar if the crossover length scales coincide, we have 
to normalize kc{£o) to use it as an independent variable 
in the range (0,1). This normalization yields kc{eo)/Tr, 
and we may see in Fig. |6] that indeed BEs and LEs fall 
very near in the first part of the Lyapunov spectrum. 
This range of approximate overlapping is limited by the 
value of n (or Eq) up to which the LVs (or the Log-BVs) 
are piecewise copies of the main LV. Note also that this 
range expands as F increases and the system becomes 
more chaotic. 

Figures [7] and [8] show that the same procedure can be 
followed also for a smaller system size {L — 128). Note 
that as the system is extensive the figures look very much 
the same as those for L = 1024. 
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FIG. 7: (color online) The colored lines with symbols indicate 
structure factors S{k) obtained for Log-BVs with different 
perturbation amplitudes for F — 6, 8, 14, 20 and L — 128. 
The values fcc(eo) are determined like in Fig. [3] 
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FIG. 8: (color online) Lyapunov spectra (solid black lines) and 
bred exponents (symbols). Size L — 128, and F = 6,7, . . . , 20. 



VI. CONCLUSIONS 

Although dynamical systems theory can characterize 
essentially all the properties of a chaotic system in terms 
of the properties of the tangent space directions (LVs) 
and their associated growth rates, the use of these lin- 
ear analysis tools in real applications is often restricted 
by practical limits. This is particularly obvious in atmo- 



spheric dynamics and weather prediction systems were 
operative models include mathematical subtleties like ad- 
hoc parametrizations of many physical processes, nondif- 
ferentiable structures, and discrete programming struc- 
tures that render the model unsuitable for linearization. 
In this context, much insight has been achieved by study- 
ing finite amplitude and truly nonlinear perturbations. 
In particular BVs have attracted much attention as a 
tool to investigate propagation of errors in both toy and 
operative weather models. However, little was known 
about the relation of finite BVs and truly infinitesimal 
LVs, apart from the obvious fact that a BV should tend 
to be coUinear with the leading LV as the BV amplitude 
£o tends to zero (independently of the norm definition 
used). 

In this paper we have uncovered a number of con- 
nections between bred and Lyapunov vectors, which are 
specially noticeable and most conveniently characterized 
when one uses the zero-norm BVs or Log-BVs. We have 
also found that these similarities appear over a spatial 
range, below some characteristic length scale, that de- 
pends on the BVs amplitude. After a Hopf-Cole trans- 
formation, the Log-BVs turn out to be a piecewise copy 
of the leading LV. This resembles what has been previ- 
ously reported for non-leading LVs [lg,[20|. Interestingly, 
the spatial structure of Log-BVs shows clearly that they 
are uncorrelated objects at long length scales over cer- 
tain characteristic length. This contrasts with LVs that 
were shown to exhibit weak correlations (decaying as 1/fc 
with the wavenumber) at long scales 16, 20]. This imme- 
diately implies that the relative heights of the plateaus 
relating Log-BVs and the LV are not constrained by the 
dynamics. In turn different initial conditions yield a di- 



versity of Log-BVs but nonetheless similar local patterns. 
Interestingly, this "regional" coincidence was observed 
before in bred vectors of a global circulation model of the 
atmosphere [25'|, but remained up to now unexplained. 
Here we have shown that this regional similarities of BVs 
of different amplitudes Sq may be understood as due to 
the overlapping with the leading LV below a certain spa- 
tial range. 

This paper has also investigated the relation between 
the growth rate of Log-BVs (the BEs) and the LEs. We 
have introduced the concept of bred exponents that de- 
scribe how nonlinear perturbations grow in time. We 
have found that BEs and LEs can be mapped onto each 
other when the crossover length scales of the correspond- 
ing vector perturbations coincide. This is only true for 
the most expanding (bred or Lyapunov) exponents, when 
the "piecewise KPZ" picture holds. The convergence of 
the BE (a finite-time LE) to the largest LE has been 
found to follow a generic scaling function, which has been 
explained by a simple argument. 

Finally, it is to be emphasized that we have used the 
Lorenz '96 model in our simulations, but very similar 
results should be expected for other dissipative systems 
with spatiotemporal chaos. 
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